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Abstract 

Wc deduce from a microscopic point of view the equation that describes how the 
state of a particle crossing a medium decoheres. We apply our results to the example 
of a particle crossing a gas, computing explicitly the Lindblad operators in terms of the 
interaction potential between the particle and a target of the medium. We interpret the 
imaginary part of the refraction index as a loss of quantum coherence, that is reflected 
in the disappearance of interference patterns in a Young experiment. 
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1 Introduction 



When a particle crosses a medium the loss of quantum coherence in its state is unavoidable. 
This means that the statistical mixture of the particle state increases as it moves on. There 
are two ways of increasing such mixture. One of them is purely classical and comes from 
a probabilistic knowledge of which the state of the medium is. The other one is specific 
of Quantum Mechanics (QM) and it requires a quantum description for the medium too. 
Here, decoherence appears because the particle gets entangled with the medium through the 
interaction with it. 

Some of these effects that the medium exerts on the particle have long been understood 
and well described by means of an effective hamiltonian, often expressed as an index of 
refraction of the medium for the particle propagation. This refraction index was introduced 
by Fermi [X in the context of neutron dispersion by matter, and is of the form: 

^ = l + 2nnl^ (1) 
k k z 

where f(k, k) is the forward scattering amplitude of the particle by a single dispersor center 
and n is the density of targets. It closely follows the calculation of the refraction index 
for light in Classical Electrodynamics when it propagates through a dielectric medium, first 
derived by H.A.Lorentz (see [2]). 

These corrections induced by the medium to the free evolution are crucial in the effects 
that matter has on neutrino oscillations |H1 E] (relevant for the MSW effect), and on kaons 

One feature of the refraction index type of correction is that the evolution of the particle 
state is still hamiltonian and, thus, its coherence is strictly preserved in time. On the contrary, 
this article is focused on the decoherence effects that are not included in a refraction index. 

A detailed study of the decoherence process would, of course, require a tracking of the 
whole system, including the environment besides the particle, which would follow a standard 
QM evolution. It is when one disregards the environment and focuses just on the particle 
state that its dynamics departs from the one described by the usual Schrodinger equation. 
It becomes non-local in time, i.e., it depends on the entire history of its evolution. Under 
certain conditions, this non-locality is very weak and decoherence can be effectively treated 
by adding new, non-hamiltonian, terms to the differential equation for the density matrix. 
(Notice that here, as decoherence is at work, the state of the system is naturally described in 
terms of a density matrix rather than a wave function). Such differentiall equation is known 
as the Lindblad equation and is of the form: 

p = -i[H eff , p]+Y^ (2A jP A] - A]A jP - pA)Aj) (2) 
j 
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The presence of an environment has two effects: it modifies in part the hamiltonian and it 
moreover generates the operators Aj which make quantum coherence of the particle itself 
decrease irreversibly. Lindblad proved equation (J2J with the assumption that the evolution 
of the density matrix in a decoherent such process fulfills quantum dynamical, completely 
positive, semigroup sort of composition law, which increases entropy [7j. This is a very 
general treatment and has been widely used to describe many different physical situations. 
Examples of these are Quantum Optics 8, 9 (where it appeared for the first time); the time 
dependence of the optical activity in chiral isomer molecules jlUj ; the emergence of a classical 
description for the macroscopic objects as well as the attempts to explain the mechanism of 
some issues of foundation in QM [111 I12j . 

Recently, it has also been the framework to address the question if decoherence could blur 
the oscillations in the experiments of neutral kaons |131 1141 115j , or the evidences of neutrino 
oscillations (^B] an d references therein). 

In this article we explore the loss of coherence for a particle as it crosses a diluted medium. 
We assume the interaction of the particle with each target is short range. For this reason 
memory effects are washed out on time-scales larger than the interaction time. We use the 
scattering S matrix to describe the interaction, since our coarsed grained time is infinite with 
respect to the intercation scale, for all purposes. A Lindblad equation thus naturally emerges. 

We discern two different sources of decoherence: foootprint and mixture, each one con- 
tributing to distinct Lindblad operators. We give explicit expressions for these operators in 
terms of the particle-target potential, their masses and the target wave functions. In partic- 
ular we compute the corrections to Fermi's equation Q in terms particle and target mass 
ratio. We also see how the unitarity of this equation is restored once all pieces of decoherence 
in the Lindblad equation are properly included. Finally, we argue that the presence of a 
medium washes out the interference pattern of a double slit Young experiment. 

2 The decoherence mechanisms: the footprint and the mix- 
ture 

We introduce, in this section, the basic ideas necessary to understand which are the micro- 
scopic mechanisms that produce decoherence in the case of a particle crossing a medium. 
We classify the sources of decoherence in two types: the footprint and the mixture. Let us 
clarify this distinction with the help of a toy model that consists of a particle and a box. 
The particle will pass through the box and will interact with it. For simplicity, suppose they 
both have two-dimensional Hilbert spaces. The orthonormal basis for the particle is {|1), |2)}, 
and for the box is {|a), After the particle crosses the box, we do a partial trace over 
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the box degrees of freedom. The two distinct mechanisms that we envisage are the following: 



The footprint. Suppose that the initial state is |l)|a). Suppose also that after they have in- 
teracted the final state is ^= (|l)|a) — |2)|6)), which is entangled. The corresponding reduced 
density matrix for the particle is the identity, which is not a pure state anymore. After the 
interaction occurs, the box "knows" what the out state for the particle is. The particle has 
left a footprint: if the box is in state \a), the particle is in |1) and if the box is in \b) , the 
particle is in |2). With the partial trace we overlook this information, that remains in the 
box as a footprint. 

The mixture. Now we start with the particle in |1), but the box in the mixed state ^|a)(a| + 
7}\b)(b\. Choose an interaction as follows: 



Notice that with this interaction and this initial state, the box does not have any "knowledge" 
about what the out state for the particle is: regardless of whether the box starts in \a) or 
\b), it remains unchanged. Yet, decoherence also appears: the reduced density matrix for the 
particle is the identity again. The source for decoherence is, in this case, the initial mixed 
state for the box . 

3 A model for the medium 

In this section we put forward the approximations that we make, as well as the procedures 
that will eventually lead to a differential equation for the particle reduced density matrix as 
it crosses the medium. 

Our medium is made of just one kind of targets. There is neither interaction between 
them nor overlapping of their wave functions. We also consider that the interaction between 
the particle and the medium is weak; therefore, we neglect terms higher than second order 
in the potential. 

Our procedure consists in dividing the medium in thin slabs of matter that will be crossed 
one by one by the particle. Each time the particle crosses a slab, we do a partial trace over the 
slab degrees of freedom and obtain a step by step equation for the particle reduced density 
matrix. Eventually, we will get a differential equation in time for it. 

The targets of each slab are in a mixed state: target j is in state \rrij) with probability 



l)|a)-|l)|a) 
1} \b) |2> |6) 



|2) \a) -► |2) \a) 
|2)|&)-|1)|&). 



(3) 
(4) 
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q m . . Since the medium is homogeneous, the states of the different targets are related simply 
by translations, and the corresponding weights are the same. 

If the initial state is pure \in) = \particle)\slab), after the particle crosses the slab the 
system is in state 

N 

\out) ~ (1 + i^T ^)\particle)\slab) (5) 
j'=i 

were T ^) is the scattering operator of the particle with target j only, and N is the number 
of targets in the slab. 

In order to justify this approximation, let us look at the case where there is just one 
particle and two targets, and consider the amplitude of the process of a particle that goes 
from x to y, and targets go from x a , Xb to y a ,Ub, in a time interval t. The terms to second order 
in the potential are those shown in Fig 1. We must integrate over the coordinates and the 
instants in which the interactions take place. Recall that the free propagator is proportional 
to the phase factor exp(im ^ 2 ^ ), m being the mass of the particle. When we integrate over 
the internal coordinates, z for instance, the oscillation of the exponential is much faster in 
picture 1.6 than in picture l.a . This is so because the size of the region where a particle 
and a target do interact is much smaller than the distance among targets. Therefore, we can 
neglect those terms of 1.6 in front of those in l.a . Ultimately, this is the same as (|5J). 




fig La fi § Lb 

Figure 1: The two order terms in V for the T matrix 

It is worth pointing out that unitary relations hold both for each 1 + iT^> and for 
1 + i~Yl!j = iT ( J ) (within the approximations above explained). Thus, terms like the crossed 
product TtWrW) are neglected by consistency. 
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4 The one step equation 



In this section we obtain the particle density matrix after it has crossed one slab. In the next 
section we will iterate this one step evolution r times. 

If the initial state for the particle is \<p n ) and for the slab is \m) = IljLi l m i)> the out- 
coming state for both together is: 

N 

|out"' m ) ~ + (j) W)H, (6) 

j=i 

The reduced density matrix for the particle is obtained by performing a partial trace over 
the slab: 

p^ = TR sta6 [|out^)(out^|] (7) 

In general, the initial state will be mixed, and will have |0") with probability p n and 11^=1 \ m j) 
with join probability n,-_i 9«v Then, an average over all possible initial states is due in (J7J). 

4.1 The footprint in the slab and the reduced density matrix for the par- 
ticle 

Let us separate from |out"' m ) the part that contains the footprint left by the particle. This 
part can be understood as the change that the particle leaves in the slab. It is thus natural 
to define it with the help of the following projector: 

P = I®|m)(m| (8) 

In order to separate the footprint term, we apply the projector (1-P) to (out"'" 1 ), which only 
keeps the part of the slab that has changed: 

(1 - P)|out"' m ) = (footprint"'™) (9) 

The part of \out n ' m ) that leaves \m) intact is: 

P|out"' m ) = \^' n ' m )\m) (10) 

One can thus write © as a sum of two orthogonal terms: 

|out n ' m ) = |</>' m )|m> + (footprint"'" 1 ), (11) 

and with it (J7J) reads: 

= (</,'"'-) (</,'"' m | + TR [(footprint"'™) (footprint"'" 1 1] (12) 
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In this way we separate the part in which the particle changes and the slab does not (the 
first term), from the one which contains any change produced in the slab (the second term). 
It is worth pointing out that the second term in Q12|) encodes the entanglement between the 
slab and the particle. From Q and (|1U|) : 

N 

= r) + i ^ K .| T o-) K)r) . (13) 

3=1 

In terms of the operators Dl?^ = (rrij\ T ^0 \mj) (which act only on the particle Hilbert 
space) the first term in (|T2|): 

^n,m ){ ^n,m l= ( 1 + i £ £>£^ |0»)<0"| (l-i^I^M (14) 

As for the footprint term, from @ and ©: 

N N 

|footprint n ' m ) = iJ2( T U) ~ ( m il T U) K> ® l{ ® N) ) 1^) II ( 15 ) 
j=i %=\ 

Each term "j" of this sum is orthogonal to any state which has target j in state \nij). For 
this reason, in the second term of (|12|) 



E TR 



1\ iV 



«g>N) 



(16) 



only the j = k terms do contribute. In order to evaluate the partial trace we introduce as a 
basis for the Hilbert space of the j particle (j!6|) then reads: 

EE (fol^W - (milT^K-Xiilmi)) |0">{0 n | (kit^^) - (m^ViX™,^)) (17) 

which in terms of the operators A^' m = T & \rrij) — (mj\T^\mj)(lj\mj) is: 

A { E j ' mj) W l ) (<t> n \A^ lj ' mi) (18) 

3,h 

Finally, we average in (|12|) over all possible initial states, for both the particle and the slab, 
with their corresponding probabilities 1 . 



{mj} 3 \ t / \ i 



(19) 



1 The unitary relations of the next section guarantee the correct normalization of p% 



out 
ed 
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where p™ ed = Yl,Pn\ ( t )n ){ c t :>n \- We call this the one step equation. Notice that the footprint 

n 

operators A^ ,m ^ are independent of the probabilities q m . . 
For the slab we obtain an equation similar to 1)19}) : 



(20) 



It is interesting to note that if p l Jl ab is the identity, and [A^ l ' n \ A^ n ^] = 0, it remains un- 
changed 2 . This may seem contradictory, because footprint effects still appear in Q19|). given 
that Ab' m ^ 0, whereas the medium has not changed. This is fixed by realizing that the 
initial state for the slab is a mixture of pure states, and the footprint is left in each one of 
them. By averaging over these pure states one erases the footprint that the particle leaves 
on the slab. Yet, the footprint effects on the particle are still there. 

4.2 Unitarity relations 

Before we finish this section, let us put forward some unitarity relations among the operators 
.A's and .D's defined above. In order to simplify the notation we omit the subindex j, and we 
write D^f and q m instead of Yli^M anc ^ II j Qmy Unitarity of S matrix reads: 

T - T f = iT^T (21) 

Let us define the operator D { ™ ] such that: T = D$ ® 1^ + D%?\ Notice that then 
(lj\DE\mj) is nothing but A l ^' mj . We can rewrite (|21|) with the help of this decomposition: 

Now let us sandwich each term of (1221) with \m). Since (m \D^\m) = 0, then: 

E A^4^ = - (i D ™ - iD\W + d\™D™) (23) 

This relation has a clear intuitive interpretation. With the approximation © for the T 
matrix, the particle cannot change the state of more than one target at a time. Thus, for a 
given initial configuration (n, {m}), the probability to change target j from state nij to state 

|| {lj \T^ \</r)\ mj ) || 2 (24) 



(22) 



2 There is similar relations between and Ab that comes from unitary relations like in 1231 . 
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The probability that any target suffers any kind of change for some initial particle state is, 
thus, written in terms of the Ae operators: 

E^EE ^ E <* n i4 (I,mj) 4 ,mj V> (25) 

n j {m 3 } Ij^rrij 

The opposite situation, in which there is no change in the slab, has a probability: 

E^E^nnKi(i+*E r(i) ) r>IlK>ii 2 (26) 

{m} n j j j 

that, written in terms of the Dm's: 

1 " E^E^I^M - D m) - V X mVmW) (27) 

n m 

These two probabilities must add up to one, which expresses the unitarity relation (J23J). 



5 The equation of motion 

In the previous section we have found an equation that relates the density matrices for the 
particle before and after it crosses a slab of matter, which we called the one-step equation. 
From it we now derive the differential equation for the evolution of this density matrix. 

5.1 From the one-step to a time equation 

Let us first justify the procedure we follow with a simple example which also applies to the 
general case. 

Suppose the particle is in a pure state and the medium is made out of static scattering 
centers. After the particle crosses the rth slab, the state for the particle is: 

\<f>(r)) = (l+iT)\<f>(r-l)) (28) 

The particle spends a time At = ^ crossing the slab, where 5 is its thickness and v is the 
mean velocity of the particle. Then, one could be naively tempted to write an equation of 
the sort: 

_ j?T| „ ~ 4 W . m (29) 



At 5 dr 5' 

But this cannot be a Schrodinger equation because T is not an hermitian operator. Demand- 
ing that (1 + iT) should be of the form exp(—iH e ffAt), the correct differential equation 
is: 

^ ) = v _ Hl + iTm (30) 
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Unfortunately, these arguments cannot be immediately generalized to the case of equation 
(|19|) . because in that case the particle gets entangled with the medium, and the one-step 
equation is much more complicated than 1)28(1 : moreover, we are aiming at an evolution 
equation of the Lindblad type, which is not a mere Schrodinger evolution. For that purpose 
lets go back to our example in (|28|) and re derive (|3Uj) in such a way that now it does generalize 
for the the one step equation ((T§|) . 

Consider that the particle crosses r slabs in a time At = r|. We can iterate equation 
(f2*5|) r times and expand (1 + iT) r : 

|#r)> = (1 + riiT) + ^^V) 2 + -MO)} 

= (1 + riiT - ^ + + O(r 2 )M0)> (31) 
Compare (|31[) with the Taylor expansion of \cp(t + At)): 

|0(t + At)) = |0(t)) + A^|0(*)) + ... (32) 

where \(p(t + At)) = \(f>(r)) and \(f>(t)) = \4>(0)). Now, identify the terms of (|31|) linear in At 
with those of (|32|) through the relation At = r-, and write: 

±\ m = V -{iT-^f + ...)m)) (33) 

which is nothing but the Taylor expansion of ln(l + iT). Thus, we retrieve (|3U|). 

We are now ready to apply this procedure to the one-step equation qi9|). To derive the 
desired relation let us Taylor expand p about the instant t, as in (|32|): 

p(t + At)~p(t) + Atp(t) + ... (34) 

On the other hand we know p after r steps by iterating equation (|19|) r times. Expanding 
the result in powers of r: 

p{r) ~ p(0) + r A [p] + ... (35) 

where A is a linear operator acting on p. By comparing equal powers of At in the two 
expressions, one finally gets: 

p{t) = - 6 K[p] (36) 

This is the equation of motion, that, as we will see in the following section, has a Lindblad 
form. 
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5.2 The equation of motion for the reduced density matrix 

Here we compute A[p] by iterating ea. l|19J) . as explained in the previous section. If we only 
keep terms linear in the potential, the iteration only gives a linear term in r: 

rE^Kp-V^), (37) 

m 

and (|36|) is an equation of Schrodinger type . It is not until we retain terms up to quadratic 
in the potential that we find any sign of decoherence and, hence, corrections to this standard 
evolution. These terms are: 

( m ) „ „•„ nt ( m ) i r>( m ) „ nt ( m ) 



El ■ n( m ) ■ nt \ m ) i 7-i( m ) r-i 

Qm [iD y M p -ip Z>y + pD M 



m 



_|_ r / r 



i) (e^mWe^ 

\ m / \ n 



t (n) 
M 



r(r ~ 1) / n (m) n (n) n tM n tW N 

m,n 

+ r £ ^4 J ' mJ V^' mj) (38) 

from which we retain only the pieces linear in r. Defining the operators A 1 ^ 71 = — Djl \ 
we can write the equation of motion (|36|) as: 

pit) = -iH P + ip rf + Mr + J E ^4^m^ } (39) 



V 

Where H = — - 
d 



E^-^Erf) 2 ) 

K {m} {m} J 



This equation is not in the Lindblad form yet. In order to do that, we use the unitarity 
relations of section (|4.2|) and rewrite it as: 



p = -i[H eff ,p] + i ^E 9mqnL^ m) [p] + ~ E (40) 

m,n lj ,rrij 

This is the final expression. We have defined the effective hamiltonian and the Lindblad 
terms as: 



H eff = ^ E 1 m M + E *n°T ) ( 41 ) 
\{m} {m} / 

r(n,m)r I _ 9 .n,m ,f n," /it n,m ,n,m A n,m »n,m 

M iPl — z M P A M A M A M P P M M 

L^ip] = 2 A%' mj) p A^' ms) - A^> m i ] A%' mj) p - p A ^h^) A (h^) ^ 
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As we can see, the hamiltonian is hermitian, and it contains the "mean" effects that the 
medium does on the particie. To first order in the potential it is nothing but the mean 
energy felt by the particle while it propagates inside the medium. 

An important remark concerning the Lindblad operators Am and Ae is in order. Mixture 
effects only depend upon how much mixed the medium is: if it is pure Am = and they are 
not present. However, footprint effects are unavoidable once the particle has modified the 
state of the medium. This always carries Ae 7^ 0, regardless of whether the state is mixed 
or not. 

5.3 The independence on 5 of the equation of motion 

Physics cannot depend on 6, the width of the elementary slab of our partition; any dependence 
on 5 in equation ()39|) is an artifact of the approximations that have been used and it will be 
automatically removed whenever the approximations hold. In order to qualify this statement, 
let us regard again, for simplicity, the case in section IS~T1 where the particle always evolves as 
a pure state: we expect that the effective hamiltonian H = | ln(l + iT^) does not depend 
on the width of the slab, while does. 

This follows from the relation between the S matrix for a slab of width 5 and the S matrix 
for a slab twice as thick: 

(1 + iT^)(l + iT&) = (1 + iT^) (43) 

from which we obtain 

^ln(l+iT( 2<5 )) = ^ln(l+iT< a >) (44) 
2d d 

This shows the independence of H on the 5 parameter. If we make some approximation in 
the computation of , the relation (|43[) will fail to hold and it will lead us to a spurious 
dependence on 5. Nevertheless, this is not substantial in the domain where the approximation 
for T is valid. 

In our case, the approximation consists of taking as T matrix the sum of the individual 
T 3 matrices of scattering of just one target j. In order for equation (|i3*)) to be satisfied: 

+ i £ T»")(l + i £ T»") = (l + i £ T*) (45) 

j€slabl j&slab2 j(Lslabl+slab2 

the crossed terms T J T l should be negligible. These are similar conditions to those imposed by 
unitarity, put forward in section (|3*|). Recall that this approximation amounts to neglecting 
re-scattering, which implies that the density of targets must be small. 
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6 The Lindblad equation for a gas 



In this section we treat the example of a gas. We expect in this case some simplifications to 
occur in the general equation (|39|) . because the gas is an homogeneus medium. In particular, 
the simplifications apply to the matrix elements (k'\H\k), that are related to the refraction 
index of the gas, as we shall see in section 7. 

6.1 The hamiltonian part 

Let us consider the hamiltonian piece from ()39[) in momentum representation: 



W«> = ^ E 1rnD^\q) + ( ~ £ q m D% 

{m} V m 




(46) 



In order to perform the computation of these sums, we divide the slab in boxes of volume 1/n 
so as to have one target inside each box (n is the density of targets). We also write \rrij,Xj) 
to label explicitly the centers Xj of the target wave functions, and similarly for the weights 
1rrij,Xj- The sum over all possible locations Xj inside boxj becomes an integral: 

(^lE^M^) = E E Im^iklD^lq) 

{m} j {rrij,Xj} 

= E E Irn, ( ndx j e i to- k >i(k\(m j) 0\T\m j ,0)\q) (47) 

We have used the fact that the matrix element (k\(mj,Xj\T\mj,Xj)\q) = e l ^ q ~ k > x i {k\(mj\T\mj)\q) , 
where it is understood that the states \rrij) are centered at the origin now. Since the weights 
q mj are the same for all targets, the sum over j, together with the integral over each box, 
becomes an integral over the slab. Then, 

(k\ ]T q mj D^\q) =n^2q m [ dxe^~ k > (k\(m\T\m)\q) (48) 

r i r J Slab 

{mj } {m} 

The integral over the slab can be split into an integral over its width (z coordinate) times 
an integral over the plane of the slab. We first perform the integral over the plane, and find 
that the result is proportional to the delta function of (q — k)», i.e., the components parallel 
to the plane: 

v r 5 

q m] D™/\q) = -(27r) 2 n Vg m / dzS^q - k)e l ^~ k ^ z {m\T(k,q)\m) (49) 
° Jo 



HI 

We compute the projection of the T-matrix in the center of mass coordinates: k cm = 
k\ + &2 and k r = ^rki — ^^2) where 1 labels the particle and 2, the target. In such 
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coordinates T reads: 



(k r \(k crn \T\q cm )\q r ) = -2ir5(k cm - q cm )5(E qr - E kr )T E (k r , q r ;m r ); (50) 

here m r is the reduced mass. The combination of these two delta functions, together with 
S»(q — k) in gives: 

S(qi-k 1 )5(q 2 -k 2 ) ,_ 1 . ■ (51) 



i fc z ,i _ k Zt 2 I 
mi ni2 I 



2m\ 



+ (5i|(gi - fci)<^ ( ki + qi - — {k 2 + q 2 ) ) <5|i(<?2 - k 2 )5 z \ k 2 -q 2 - 2qi + — (k 2 + q 2 ) . , 

The delta functions in the two terms provide the energy momentum conservation laws corre- 
sponding to one-dimensional elastic scattering. The first term is diagonal in both particle and 
target momenta; it corresponds to the process where each velocity remains unchanged and 
gives the dominant contribution to (|49|) . because the imaginary argument of the exponential 
vanishes. The second term corresponds to the solution where the relative velocity flips its 
sign, and gives a negligible contribution if the particle wave packet is entirely contained inside 
the slab of width 5: therefore k z 5 ^> 1 and the oscillations of the exponential factor kill the 
integral over z. Then (|49j) reads 3 : 

^{h\ ^2 QmjD^/lqi) = -(2TT) 3 n^2q m {m\- — j— T E (k r , k r ; m r )\m)S(ki - qi) (52) 

{rrij} rn I m2 k z X I 

The hamiltonian (|4*H)) contains, also, the square term (%2 m q mj D™j ) 2 . We discard this 
piece, because it is quadratic in n, which should be small by consistency of the approximation 
(see (|45p ). Finally, we get for the H matrix elements: 

(ki\H\qi) = (27r) 3 nV )q m (m\ — r — T E (k r ,k r ;m r )\m)d(ki - q{) (53) 

£ — ' I -i mi >^z,2 I 

m I 1 ~ m 2 fc z ,i I 

We assume that the tipical velocity of a target is much smaller than the particle velocity, and 
therefore the denominator |1 — lULtiAl ^ i regardless of whether — is small or not. 

I m 2 fc z ,i I » o m2 

This is our final expression for the hamiltonian of the particle. As we can check, there is 
no-dependence on the width 5 of the slab left over in 1)53(1 . Further simplifications occur in 
the limits of both mi m 2 and m\ 3> m 2 , that we now consider. 

The limit of a heavy target. 

Let us expand (|52*1) in powers of ^ <C 1. The zeroth order term is: 

{ki\B\ qi ) = (27r) 3 nS(k 1 - q^TEihM^ + of^) , (54) 



3 Recall that = v, since we assume that particle state is peaket around velocity v pointing in z direction 
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which is independent on the state of the targets. 

There are two sources of corrections to this: the denominator of (|52|). and the T E matrix 
elements through their dependence on k r and m r . In a gas, the target momentum expectation 
value vanishes, i.e., qm(m\k2\m) = 0, and the corrections coming from the denominator 
cancel. The first — correction, thus, reads: 

- S(ki - qi)(2n) 3 n — ki—-T E (ki,ki;m r ), (55) 

7TT-2 OK\ 

where we have not expanded m r in T E . Notice that these corrections start quadratic in the 
potential because the linear term in V does not depend on k r at all. The final expression in 
this limit reads: 

(h\H\ qi ) = 5(h - (? i)(2vr) 3 n ( 1 - ^(k 1 ^-)) T E (k u h;m r ) + O ( (—) 2 ) (56) 

V m 2 oki J \ m 2 J 

The corrections to the infinitely heavy target limit are 0(^, V 2 ). 
The limit of a heavy particle. 

Let us now consider the opposite limit, and expand ()52|) in powers of ^ <C 1. The zeroth 
order term is: 

(fei|IT|gi) = (2vr) 3 n V q m (m\T E (^h - k 2 , - k 2 ;m 2 ) \m)d(h - ft) + O (—) (57) 

L — ' \ mi mi / \ mi / 

^■ki = m 2 v\ should not be treated as a small quantity as compared to k 2 if we assume that 
the particle velocity v\ is larger than the velocities of the targets in the gas. 

The first order corrections contain a term coming from the denominator of (|52|) which 
does not cancel as it did before. Finally, the hamiltonian (k\\H\q\) reads: 

(2vr) 3 n V q m (m\(l + + k 2 JL))T E (K = - kK; m r )\m)S(h -q x ) + (( — ? 

z — ' mi kr oK mi V mi 

m x 

where we have substitute the quotient k z ^ 2 jk z ^\ by k\k 2 jk\: they both coincide in the limit 
of a particle state pecked in the z direction. 

6.2 The piece of decoherence 

Apart from H, the evolution equation (|39j) involves further pieces that contain the Am and 
A E operators, which are the ones that induce decoherence. Let us first consider the terms 
containing Am (see section ()5.2|) for details ): 



It — V X - ^ A n,m t]n,m v V - ^ f-.m rit m V r\m \ ^ nt n 

M = 2 J <7m<?n^ M pA^ = - 2^ q m D M pD M -gl^ <1™ D mP 2^ ^n-DV 
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Making explicit the summation over all targets in (|59|) we find: 

v 
~6 



A MP Al = ^ E E " 9 2 mj , Xj )DZ j ' Xj pD^ (60) 



where we have split label (nij) into (nij,Xj) as before, to explicitly display the center of target 
j wave function. We proceed as before and convert the sum over Xj into an integral (note 
that the squared probability in (|6U|) vanishes). We get: 

A M pA j M = E 9m / d 3 xD^ x pD^ x (61) 

m J slab 

The term coming from the Ae operator will have a similar expression, yielding: 
*bpA b = - s E ^4 J ' mJ V^' mj) =njE^E / ^4™^'™'* (62) 

m^lj^mj m l^m J slab 

The matrix elements (&', /c) of these two terms together add up to: 

(k'\A MP A ] M + AjspAklA;) = n| E 9m J dqdq' {2v) 2 5\\ (k' + q - k - q') x (63) 



I, in 

s 



x / dze i ^- k ^-^- q ^ z {k'\(l\T\m)\q')p(q',q) (q\(m\T*\l)\k) 
Jo 

The T-matrix elements involved in the last expression, to first order in the potential, and in 
momentum representation, are: 



(fci|<i|r|m>|^) = -2vry J d 3 k' 2 d 3 q' 2 AHk' 2 )5(k' cm - q' cm )8{E q , - E k ,)V(q r - k r )A m (q 2 ) 
(*i|<l|Tt|m)|gi) = -2tt / I d 3 k 2 d 3 q 2 A l {k 2 )5{k cm - q cm )5(E qr - E kr )V*(q r - k r )A m (q 2 )(M) 



e i((fci-*,)-(gi- 



where A[, A m are the target wave functions. A resolution of the identity ^ A^(k' 2 )Ai(k 2 
5{k 2 — k 2 ) appears in ((51)) once expressions ((51)1 are used. 

(k'\A MP A M + A EP A E \k) = n^(27r) 4 E<7m J d 3 k 2 d 3 qxd 3 q 2 d 3 k' 2 d 3 q[d 3 q 2 J dz 

l,m 

x5(k 2 - k' 2 )8\\ (k[ +qi-h- q[) 5(k cm - q cm )8(k' cm - q cm )5{E q > r - E K )5{E q 
xA m (q 2 )A* m (q' 2 )V*(q r - k r )V{q' r - k' r ) p(q[, Ql ) 

This is our final expression for Lindblad piece of the equation of motion (|39|) . 

Further simplifications occur, again, in the limits of heavy target and heavy particle. The 
essential difference between the two limits enters through the energy delta functions. 

In both cases, the following relation is useful to write a perturbative expansion in the 
ratio of masses: 

5(E qr - E kr ) = -L ( 1 - lr—-) (66) 

2m \E qr - E kr - le E qr - E kr + it ) 
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We now study the two different limits separately. 



The limit of a heavy target. 



In such limit we can expand the terms in (|66|) as a geometric series: 



^—T2=- (67) 



E qr -E kr ^ie q\ - k\ =fie ^ \2 m t q\-k\^-ie 

with the caution that with this expansion the series thus obtained is asymptotic. Here we 
have used the conservation of the center of mass momentum 5{k cm — q cm ) present in (|65|) . 
By taking higher derivatives of (|66|) we have 



1 

27d 



(-iyd s 6(E qi -E ki : 

si dE° qi 



(68) 



(E qi - E kl - ie) s+1 (E qi - E kl + ie) s+l 
that allows us to write (|66j) as: 

S(E qr - E kr ) = ± ^/ dS6{E « E ; E ^ (fe + - (q 2 + kl )r (69) 
s=0 31 

suitable for an expansion in inverse powers of the target mass. Let us focus on the zeroth 
order term of this expansion. The set of delta functions in (|65|) is now: 

5{k 2 - k' 2 )8\\ (k[ +qi-h- q[) 8{k cm - q C m)S(k' cm - q CTn )5(E qli - E kli )8{E qi - E kl ) (70) 

At this point, we replace such string of deltas by the following expression 

■5(k 2 - k' 2 )5 (k[ +qi-h- q[) 5(k cm - q cm )5{k' cm - q' cm )5(E qi - E kl ) (71) 



\ v z,l 



Let us justify this substitution. In first place, the two terms have the same trace in the limit 
where the state for the particle is peaked around a velocity v and the size of the wave packet 
is smaller than 5, the width of the slab. Also, in such limit (where p(q, q') is different from 
zero only if q ~ q'), the regions of momenta q,q',k,k' for which the arguments of the delta 
functions vanish are the same in both expressions. In this limit of q « q', from (|7U|) we find 
in principle two allowed regions: 

(k' - fe)|| =(q- gOu w ; E k , = E q , « E q = E k (k ~ k') or (k[ - k\\ , k z ~ -k'^j (72) 

but the only one that survives upon z integration in (|65|). is the first one, with k ~ k' . This 
is the same region that makes the arguments of the delta functions in (|71|) vanish too. 
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The zeroth order term in \/m 2 of l|65|) is, thus: 

(k'\A M pA M + A EP A E \k) = (73) 
n(2vr) 4 y d^dg^ (*£ + 9l - A* - <£) 5(E qi - E kl )V*{ qi - h)V(q[ - k[) p^, q x ) 

Notice that this is totally independent on the specific form of the target states and preserves 
the mean value of the particle energy. 



The limit of a heavy particle. 



Let us now consider the opposite limit of a heavy particle. The terms in ()66|) can also be 
expanded as a geometric series, now in powers of The relative energy E qr is: 

E gr = + ^ + ^q 2 ? (74) 

2m r m t mt 

As in the hamiltonian piece, we cannot treat ^qx as small quantity in front of q 2 , since we 
assume that the particle velocity V\ is larger than the velocities of the targets in the gas. 
We rewrite E qr as: 

E Vr = ^—(-Q2 + m rVi + — (Agi + q 2 )f (75) 
2m r mt 

where Aq\ = qi — mivi. This brings our expressions to a form that can be easily expanded 
in powers of ^(Aqi + q 2 ), the fluctuations of the particle and target momenta about their 
approximate mean values. (|66|) reads: 

1 2m r ^ ( 2m 2 (k 2 - tfeXAgi + g 2 ) 

s=0 



y . 

E qr - E kr =F it {m r vi - q 2 ) 2 - (m r vi - k 2 ) 2 =F «e \ m t (m r vi - q 2 ) 2 - (m r vi - k 2 ) 2 =p «e 



that is an expansion in inverse powers of the particle mass. 

To zeroth order, the delta functions that appear in ([65(1 are now: 

<5(fc 2 — ^2)^11(^1 + <7i — ^1 — q'i)3(kcm — <lcm)5{k' cm — q cm )5(E m2V1 _ q2 — E m2Vl _ q t i )5(E m2V1 _fa — E m2Vl _ q2 ) 

At this point we use the fact that the velocity of the particle is much bigger than the target 
velocity, and the last expression reads: 

-— -5{k 2 - k' 2 )5(q 2 - q' 2 )8{ k l + <?i ~ h ~ + k 2 - qi - q 2 )S(E m2Vl ^ k2 - E m2Vl _ q2 )(78) 

\v z ,l I 

Finally, equation (|65f) is, to zeroth order: 

(k'\A M pA\ I+ A EP A E \k) = n(2vr) 4 ^ q m f d 3 k 2 d 3 q 2 d 3 q[d 3 q x (7 
5{k[ + q X -k x - q[)5(k x +k 2 -qi- q 2 )5(E m2Vl _ k2 - E m2Vl _ q2 )p(q' 1 , q x )\A m (q 2 )\ 2 \V* (k 2 - q 2 
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7 The index of refraction 



The index of refraction is usually defined as the phase shift that the medium induces on a 
plane wave as the particle propagates through the medium. In an evolution where coherence 
is not preserved, this concept, as it stands, does not hold anymore. This is the case of the 
Lindblad evolution. Yet, a generalization is still possible for the part of p that preserves 
coherence. Let us first identify such part regarding the way a Lindblad equation increases 
the mixing of the density matrix. A Lindblad evolution is of the form: 

p = -iHp + ipH ] + 2ApA\ (80) 

where H is H e ff — iA^A. After an infinitesimal time interval, the density matrix evolves to: 

p{t + St) = p(t) + iSt(Hp - pH^) + 5t(2ApA^) (81) 

The first two terms together add up to a density matrix which has trace less than one (since 
H is non-hermitian) and coherence is preserved on it (since it is an hamiltonian evolution). 
There is a further term, the last one in (|81|). which is also a density matrix on its own. Thus, 
we see that under such a time evolution a density matrix can always be split as a part in 
which coherence is preserved, p coh , plus the rest, p mix , which is also a density matrix (neither 
one is normalized here). This mechanism reflects the fact that the average of two density 
matrices is, in general, a density matrix with more mixing. This suggests a definition of p coh 
so that it satisfies the hamiltonian equation: 

d 



dt P 



~ coh - iHp coh + ip coh H^- (82) 



then, for p r ' 



^_ p mix = _ iHp mix + ip mix H ] + 2 Ap mix A^ + 2Ap coh A\ (83) 

which is the same as (|80|) plus a term that depends on p coh . This non-homogeneous term is the 
responsible for the increasing of the trace of p mM! , in order to maintain properly normalized 
the total density matrix p. Equation (|83|) also ensures that p mtx is a positive matrix, since 
its time derivative is a sum of positive matrices. 

It is natural, thus, to write p in equation (|39j) as p coh + p mix , with p coh and p mix satisfying 
the equations: 



— n coh 

dt p 



HpCoh _ p coH H ^ ffE Wj qmD ^ _ qmD ^)f ( 84 ) 

" w / 

r\ -1-1 

® n mix/j.\ • f tt mix 1 , ^ \ r j (n,m) r mill , ^ \ j (h > m j) r mixi 



u n mixi.\ _ ■ \ tt mixl i x v _ n T{n,m)t mixi , x v \ " n j 

—p (t) _ -i[H eff ,p J + -- ^q m q n L M [p I + ^t^M 



711,71 



--S^ a n A n,m n coh A 




IP 



lv Sr- n „ A n,m coh A ]n,m v . (!j,mj) n coh A(h > m j) /oe\ 
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We naturally take p coh (0) = p(0) and p miX (0) = as initial conditions. 

In spite of the fact that p decoheres, an index of refraction can be defined from p coh , 
since it still evolves with an hamiltonian, which has already been computed for the case of 
an homogeneous gas in terms of the forward scattering amplitude. In order to retrieve the 
refraction index we should add the free hamiltonian, and get: 



k> I (k\H\ 



k y k 2 /2m\ 
In the limit of heavy targets, ^ < 1, we obtain: 



(86) 



k' . , , „ mi (l-^kJL))f(k,k,m r ) / mi 



« i + 2nn^- + o (^f (87) 

k m r fc z \ m>2 J 

This is the generalization of the result obtained by Fermi in the case of scattering centers . 
Notice that the imaginary part of (|87|) has a clear interpretation: it is the depletion of the 
amplitude modulus due to a loss of coherence by dispersion. 

The opposite situation of heavy particle, when 2^ 3> 1, has also been computed in the 
previous section, with the resulting index of refraction: 

k' mi C 

— « l + 2vrn— 

k m r k z 

where the function C is: 

C = V<? m (m| (l + —(hh/kf + k 2 ^-)) f(K = —ki - k 2 , K;m r )\m) 
^-^ \ mi oK J mi 

8 Decoherence effects and interference patterns 

The decomposition found in the previous section allows us to see how interference patterns 
are destroyed by the decoherence induced through the interaction with the medium. Such 
patterns appear by the interference of two states, \<f>) and that travel following different 
paths and rejoin later. Thus, a superposition of \<j>) and \ip) must retain enough coherence 
so that these interference patterns can be observed. For the rest of this section we assume 
that the initial state |<E>) is the sum \<jj) + We thus start out with p coh (0) = \&)(&\ and 
p mix (0) = 0. 

The crossed terms of the p matrix are the responsible for the interferences: we argue 

that in the case when \<f>) and are well separated and localized they decay exponentially 
in p(t). The proof is in two steps. We first recall that the whole p coh decays in time, 
straightforwardly from its construction, in particular its crossed terms. Then we see that 
these crossed terms never feed-back into p through p mix via eq. (|83|) if the two states are 
separated enough and localized. 
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The proof goes by integrating the differential equation (|83|) for p mtx for an infinitesimal 
time interval 5t: 

p ™* m = 5t l -- 5 Y, qmqn Airp coh mtr + «| E Q m A j,mj) p coh (v)4 li,mj) m 

m,n L.rrii 



All crossed terms \4>)(ip\ appearing on the r.h.s. of this equation vanish. The crossed terms 

(l- m ) 

in the second piece are directly zero, because A E J ' 3 is the amplitude for target j to jump 
from \rrij) to by interacting with the particle: when target j is close to the state \<p) it 
is far from the state then A^' m = 0, and vice-versa. As for the crossed terms from 
the first piece, they are proportional to: 



E lift* IK - d%) E(< 1 - d m) (9i) 

{m,n} \ I I' J \ j j V i J 

and since states \4>) and \~ijj) are far away, targets that contribute in the summation over the 
j index are different from those that contribute in the summation over i index. This allows 
us to write (|91|) as: 

E E <l m] q nj (D^ - D n A \4>)(M I E E ImSnMr - D^)\ =0(92) 
v {j close to cf>} m i> n i J \{i close to V}" 11 '™ 1 / 

Then, after an infinitesimal time 5t, crossed terms stay out from p mtx {t). This completes the 
proof. 

We apply this result to the study of the interference patterns in a two slit Young experi- 
ment. We shall study separately the probability distribution of hits on the screen that comes 
out from p mix ^ and from p coh , each one being a density matrix on its own (except for an 
unessential normalization factor). Finally, the observed pattern is retrieved as the sum of 
both. 

The part of p coh is easy: It is essentially the same as if there were no medium, with 
the wave number corrected k — > k' . This changes the oscillation wave-length to % 
and introduces a global damping factor of e - 2Im ( k ') L which diminishes the amplitude on the 
screen by this amount. 

Unlike p coh , the part of p mix is far more involved to obtain, for it contains all the de- 
coherence effects and it is the one which eventually washes out the interferences when the 
propagation of the particle between the slits and the screen is through a medium, rather 
than the vacuum. However, it is easy to see that it produces an approximately constant 
background in the region close to the central peaks. 

The argument is indirect and uses the result of an auxiliary problem, that has as initial 
state p(0) = gl^K^I + ^IV'XV'I; f° r which the solution can be qualitatively described. For 
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such a totally uncoherent state there are no oscillations and the intensity is the sum of the 
intensities produced with one of the slits closed. In the presence of matter the curve is rather 
broad and flat, because each target disperses the particle in all directions, thus increasing 
the odds that it hits the screen away from the center. Such a broad shape is only possible if 
the mixing piece itself is even broader, since p coh only reduces its size by a depletion factor 

e -2Im{k')L 

This is already the end of the argument if one realizes that this p m%x has to be the same 
either with the initial auxiliary state p(0) or with \<p) + because - and this is crucial - 
the crossed terms of the initial state do not intervene at all in the subsequent p mtx 

evolution, as previously shown. 

We thus see that the mixing part adds nothing but a constant background around the 
central region, that dies off very slowly as one moves away. 

Summarizing, the observed interference pattern consists of a constant background, which 
entails the decoherence due to the medium, superposed to the oscillations of p coh , also damped 
by the medium (see Fig. 2). The relative strength of the two contributions is dictated 
by unitarity, and one can roughly estimate the ratio of the oscillation amplitude over the 
background as: 

2 e -2Im(k')L 

(93) 

i.e., the smaller this number may get, the more invisible the interference fringes become. 

Finally, we would like to stress the observation that from the size of the tiny oscillations of 
the interference pattern in the central region, one can fit e ~ 2Im ( k ') . Precise fits in experiments 
of this kind thus provide a direct measurement of Im(k'), related to Im(H) = —A^A in eq. 
(|SU|) ; i.e., the size of the Lindblad coefficients for the particle crossing a medium. 

9 Conclusions 

In this article we have addressed a microscopic derivation of the Lindblad equation that 
describes the loss of quantum coherence of a particle that crosses a medium with which 
it interacts. This is done for a dilute medium made out of targets that interact with a 
particle through a short range potential. In the case that the medium is a gas, we keep in 
the evolution equation terms up to second order in the potential. We also use the impulse 
approximation, i.e., we approximate the scattering amplitude for the particle with the gas by 
the sum of the scattering amplitudes with each single target. We argue how the imaginary 
part of the refraction index given by Fermi is related to the Lindblad terms and is, therefore, 
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Figure 2: The interference pattern wen a gas is present: 

due to decoherence. Finally, we show that when we apply our results to a double slit Young 
experiment, the characteristic interference patterns disappear due to the presence of the 
medium, a distinctive sign of decoherence. 

We stress how the corrections due to the presence of matter split into different terms ac- 
cording to the two different mechanisms that produce decoherence: the footprint mechanism, 
related to the entanglement with the medium; and the mixture mechanism, related to the 
probabilistic knowledge of the medium state. 
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